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Transport  through  ionic  conducting  membranes  is  examined.  An  equation  describing  the  chemical  poten¬ 
tial,  fis,  of  electrically  neutral  species,  s,  in  the  membrane  is  derived  in  terms  of  ionic  and  electronic 
currents,  and  ionic  and  electronic  transport  resistances.  It  is  shown  that  the  fis  in  the  membrane  need 
not  be  mathematically  bounded  by  the  values  at  the  two  electrodes  (reservoirs)  if  the  ionic  and  the 
electronic  currents  through  the  membrane  are  in  the  same  direction.  Conditions  could  develop  under 
which  the  fis  in  the  membrane  may  exceed  the  thermodynamic  stability  of  the  membrane  even  when 
exposed  to  stable  conditions  at  the  two  electrodes.  It  is  shown  that  during  charging,  chemical  potential  of 
lithium,  fiu,  in  the  electrolyte  of  a  lithium-ion  battery  may  exceed  that  corresponding  to  pure  lithium  thus 
causing  lithium  precipitation  and/or  reaction  with  the  electrolyte.  It  is  also  shown  that  in  a  lithium  ion 
battery  pack  containing  several  cells,  degradation  may  occur  during  discharge  due  to  cell  imbalance.  In 
unbalanced  cells,  the  SEI  layer  may  form  at  both  the  anode/electrolyte  and  the  cathode/electrolyte  inter¬ 
faces.  A  bi-layer  separator  comprising  an  electronic  conductor  and  an  electronic  insulator  is  proposed  for 
improved  stability  of  lithium  batteries. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

It  is  well  known  that  many  electrochemical  devices  such  as  bat¬ 
teries,  fuel  cells,  and  electrolyzers  degrade  over  time  and  under 
various  operating  conditions.  In  batteries,  for  example,  it  is  known 
that  in  series-connected  cells,  degradation  often  occurs  if  one  or 
more  cells  exhibit  different  characteristics  than  the  rest  of  the  cells 
[1,2].  This  phenomenon  in  batteries  is  attributed  to  ‘cell  imbalance’, 
and  great  care  is  necessary  to  minimize  failures.  Cell  imbalance  has 
been  extensively  investigated,  and  various  engineering  solutions 
have  been  developed  [3].  Many  are  based  on  the  use  of  balanc¬ 
ing  circuits.  Also,  batteries  are  more  prone  to  degradation  during 
charging.  In  lithium  ion  batteries,  the  design  of  separator,  the  choice 
of  electrolyte,  additives  to  electrolytes,  the  incorporation  of  sen¬ 
sors,  etc.  are  some  of  the  approaches  used  to  increase  safety  [4,5]. 
It  is  known  that  fuel  cells  (especially  stacks)  also  undergo  degra¬ 
dation  during  operation.  There  is  considerable  literature  on  solid 
oxide  electrolyzers,  which  too  are  known  to  degrade  under  certain 
operating  conditions  [6-8].  A  review  of  literature  shows  that  many 
electrochemical  devices  degrade  under  various  operating  condi¬ 
tions.  While  there  are  several  reasons  for  degradation,  and  many 
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depend  upon  the  specifics  of  a  particular  system,  the  observation 
that  many  electrochemical  devices  degrade  suggests  that  there 
may  be  a  common  underlying  reason  (or  reasons)  which  is  (are) 
applicable  to  all  such  electrochemical  systems.  The  intent  of  this 
manuscript  is  to  propose  a  mechanism  of  degradation  based  on  fun¬ 
damentals  of  linear  non-equilibrium  thermodynamics  which  may 
be  applicable  to  virtually  all  active  (current  carrying)  electrochem¬ 
ical  devices. 

An  active  electrochemical  device  such  as  a  fuel  cell,  a  battery,  and 
an  electrolyzer  involves  transport  of  ionic  and  electronic  species 
under  thermodynamic  potential  gradients  (forces).  Much  of  the  lit¬ 
erature  on  transport  is  based  on  expressing  fluxes  as  functions  of 
forces.  The  simplest  case  is  that  of  transport  of  a  species  in  response 
to  a  force  acting  on  it.  Many  phenomena,  however,  involve  coupling 
of  two  or  more  processes  [9,10].  For  instance,  transport  of  two  or 
more  species  may  be  coupled  [9,1 0].  In  many  such  cases,  fluxes  can 
be  written  as  linear  functions  of  forces  and  transport  coefficients. 
The  fundamental  basis  of  linear  non  equilibrium  thermodynamics 
applied  to  coupled  processes  is  the  classic  work  by  Onsager,  accord¬ 
ing  to  which  the  generalized  flux  j;  of  species  i  is  linearly  related  to 
the  generalized  force  Xk  on  /<,  by  the  equations  [9-11] 


ji=j2LikXk  (i) 

k 
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Nomenclature 

di 

thermodynamic  activity  of  i 

e 

electronic  charge  (C) 

£n 

Nernst  potential  (V) 

Ea 

applied  voltage  (V) 

EceU 

average  cell  voltage  (V) 

V 

voltage  per  cell  (V) 

X 

position  (m) 

h 

ionic  current  density  through  the  membrane 

(Am-2) 

h 

electronic  current  density  through  the  membrane 

(Am-2) 

Lik 

Onsager  transport  coefficients 

Rik 

Onsager  transport  resistances 

N 

number  of  series  connected  cells  in  a  battery  pack 

ji 

thermodynamic  flux  of  species  i 

Xk 

thermodynamic  force  on  k 

Ri 

membrane  ionic  area  specific  resistance  (including 
interfaces)  (£2m2) 

Re 

membrane  electronic  area  specific  resistance 
(including  interfaces)  (£2  m2) 

r 

position  vector 

s 

neutral  species 

T 

temperature  (K) 

n(x) 

membrane  ionic  area  specific  resistance  between  a 
point  in  the  left  reservoir  (electrode)  and  position  x 
in  the  membrane  (£2  m2) 

re(x) 

membrane  electronic  area  specific  resistance 
between  a  point  in  the  left  reservoir  (electrode)  and 
position  x  in  the  membrane  (£2  cm2) 

m\,m2 

masses  of  spheres  (kg) 

h/2 

radii  of  spheres  (m) 

V\V2 

velocities  of  spheres  (m  s-1 ) 

g 

gravitational  acceleration  (ms-2) 

T 

tension  in  the  string  (N) 

diS/dt 

entropy  generation  rate  (J  K-1  s-1 ) 

Zi 

valence  of  i 

<5*o 

interface  thickness  atx0  (m) 

** 

interface  thickness  at  Xf  (m) 

Mi 

chemical  potential  of  i  (J  species-1 ) 

Ai 

electrochemical  potential  of  i  (J  species-1 ) 

Mu 

chemical  potential  of  lithium  (J  species-1 ) 

standard  state  chemical  potential  of  pure  lithium 
(J  species-1) 

0 

electrostatic  potential  (V) 

<P 

electric  potential  (V) 

&i 

ionic  conductivity  (S  m-1 ) 

<?e 

electronic  conductivity  (S  m-1 ) 

Tf] 

viscosity  of  the  fluid 

V 

dissipation  function 

where  Lik  are  the  transport  coefficients.  Onsager  reciprocal  rela¬ 
tions  show  that 

Lik  —  Lki  (2) 

The  coefficients  I2  ^  k  are  the  coupling  coefficients  (non-diagonal 
terms  in  the  [Lik]  matrix).  For  a  two  species  transport,  the  equations 
are 

h  =LUX1+LUX2  (3) 

h  =L2\X\  +  L22X2  (4) 


If  the  forces  are  independent  of  each  other  and  the  fluxes  are 
independent  of  each  other,  it  is  known  that  [10] 

LuL22-L2u>0  (5) 

This  inequality  is  called  positive  definite.  The  two  limiting  cases 
of  Eq.  (5)  are  as  follows:  (a)  L\\L22  -L^2  =  0  (and  thus  actually  not 
positive  definite).  This  case  corresponds  to  full  coupling,  (b)  l\2  =  0, 
which  corresponds  to  no  coupling. 

In  most  studies,  certain  experimental  conditions  are  imposed 
that  create  forces,  Xks,  which  drive  the  fluxes  J/S.  The  forces  Xks  are 
not  easily  measurable,  but  the  fluxes,  J2s,  are  usually  measurable. 
Thus,  in  the  often  used  formulation  of  flux-force  relationships  (Eq. 
(1)),  the  easily  measurable  parameters  are  on  the  left  hand  side  of 
the  equations  while  those  not  so  easy  to  measure  are  on  the  right 
hand  side.  Onsager  gave  the  equations  in  both  forms;  fluxes  written 
as  functions  of  forces  and  forces  written  as  functions  of  fluxes.  The 
latter  are  given  by 

*i  =  (6) 

k 

in  which  the  measurable  parameters  are  on  the  right  hand  side 
(independent  variables)  and  difficult  to  measure  parameters  are  on 
the  left  hand  side  (dependent  variables).  The  reciprocity  relations 


give 

Rik  =  Rki  (7) 

For  a  two  species  transport,  the  corresponding  inverted  equa¬ 
tions  written  with  fluxes  as  independent  variables  are 

Xi  =Ruh+Ri2l2  (8) 

%2  =  K21J1  +  R22J2  (9) 

and 

*11*22 -*?2>0  (10) 


The  relationships  among  the  various  Onsager  coefficients  (Ly-s 
and  RyS)  depend  upon  whether  all  of  the  forces  are  independent, 
all  of  the  fluxes  are  independent,  or  whether  some  of  the  forces 
and/or  fluxes  are  dependent.  Hasse  [10]  has  discussed  some  specific 
examples  for  two  and  three  species  systems. 

For  matter  transport,  the  forces  can  be  related  to  gradients  in 
chemical  potentials.  Thus,  Xk  is  related  to  the  chemical  potential  of 
k,/A,k  by  [10] 

=  (11) 

and  thus 

Ji  =  -5>V/,fc  (12) 

k 

The  preceding  equations  for  two  species  transport  thus  are 


J\  =  -fnV/zi  -L12V/x2 

(13) 

h  =  —L2\  V -L22V/X2 

(14) 

and 

Vmi  =  -RnJi  ~  R\2J2 

(15) 

VM2  =  -^21/1  -  R22J2 

(16) 

Both  forms  of  these  equations  have  been  used  in  the  litera¬ 
ture  and  discussed  in  numerous  textbooks.  However,  the  next  step 
which  relates  the  chemical  potentials  to  fluxes  does  not  appear  to 
have  been  discussed  in  the  literature,  at  least  not  from  the  stand¬ 
point  of  stability  [12].  The  obvious  next  step  is  the  integration  of 
forces  given  as  functions  of  fluxes  to  obtain  the  chemical  poten¬ 
tials.  While  mathematically  a  trivial  step,  it  is  the  intent  of  this 
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manuscript  to  propose  that  this  next  step  is  important  from  the 
standpoint  of  stability  of  the  materials/systems  ‘subjected  to  trans¬ 
port’.  The  chemical  potential  of  species  i  as  a  function  of  fluxes  and 
spatial  coordinates  is  then  given  by 

Mi(?)  -  Mi(?o)  =  r0)  =  -  f  y^Rik]k  ■  dr  (17) 

Jf„  k 

where  r  and  r0  are  position  vectors  from  some  arbitrarily  cho¬ 
sen  origin.  Thus,  Eq.  (17)  gives  differences  in  chemical  potentials 
of  species  i  between  two  positions  (r  and  r0).  Eq.  (17)  states  that 
if  the  fluxes  of  all  species  are  known  as  functions  of  position,  if 
all  (inverse)  transport  coefficients  are  known  as  functions  of  posi¬ 
tion,  and  if  the  chemical  potentials  of  all  species  are  known  at 
one  position  (r0),  then  the  chemical  potentials  of  all  species  can 
be  determined  at  any  other  position  (r)  over  which  the  transport 
parameters  (Ly  s  or  Ry  s)  are  known.  In  this  manuscript  we  will  only 
discuss  the  steady  state,  thus  no  time  dependence  is  involved.  The 
approach,  however,  can  in  principle  be  extended  to  include  the  time 
dependence  (transient  cases). 

If  we  know  all  of  the  fluxes,  no  matter  what  their  magnitudes, 
we  know  nothing  about  whether  or  not  a  system  is  thermodynami¬ 
cally  stable,  since  transport  equations  involve  gradients  in  chemical 
potentials— not  chemical  potentials.  However,  in  an  integrated 
form,  Eq.  (17)  provides  magnitudes  of  the  chemical  potentials  of 
various  species,  which  is  the  information  needed  to  determine 
stability  under  transport.  A  system  (or  a  material)  is  stable  over 
some  given  chemical  potential  range  of  a  species  (or  given  chem¬ 
ical  potential  ranges  of  more  than  one  species)— but  outside  of 
the  range,  the  system  (material)  may  not  be  thermodynamically 
stable. 

The  objective  here  is  to  determine  what  may  be  the  ‘externally’ 
imposed  conditions  of  fluxes  of  the  various  species  on  a  system  (a 
separation  membrane,  a  fuel  cell  membrane,  a  battery  membrane, 
an  electrolyzer  membrane,  etc.)  which  would  force  the  chemi¬ 
cal  potentials  of  one  or  more  (electrically  neutral)  species  in  the 
membrane  go  outside  the  thermodynamic  stability  range  of  the 
membrane.  Fig.  1  shows  a  schematic  and  the  significance.  With  ref¬ 
erence  to  Fig.  1 ,  our  focus  is  the  membrane  between  x0  and  Xf ,  for  the 
time  being  without  discussing  the  details  of  interfaces  separating 
the  membrane  and  the  reservoirs  (electrodes).  There  are  two  tran¬ 
sition  regions  at  the  two  interfaces  over  which  properties  change 
sharply  over  a  very  small  distance  (perhaps  a  few  nm  or  even  a 
fraction  of  a  nm).  In  Fig.  1,  the  insets  show  enlarged  schematics  of 
interface  regions  of  thicknesses  8Xo  and  8Xf .  Thus  on  a  microscopic  or 
even  on  a  submicroscopic  scale,  properties  change  abruptly  across 
the  interfaces. 

The  problem  examined  here  is  different  from  those  given  in 
books  on  irreversible  thermodynamics  (e.g.  Hasse  [10],  Denbigh 
[11]).  For  example,  Hasse  [10]  has  examined  transport  across  an 
interface  (a  thin  membrane)  separating  two  phases  (reservoirs). 
The  focus  of  most  such  studies  is  on  transport  between  two  reser¬ 
voirs  through  a  membrane  but  not  on  the  ‘state’  of  the  membrane. 
In  those  studies  which  have  focused  on  the  membrane  (and  there 
are  numerous),  the  emphasis  has  been  on  transport  but  not  on 
membrane  thermodynamics.  The  present  focus  is  on  membrane 
thermodynamics  (that  is  the  ‘membrane’  separating  the  two  reser¬ 
voirs  or  electrodes  is  the  ‘system’  of  interest).  In  studies  reported 
to  date,  it  has  been  assumed  (explicitly  or  tacitly)  that  chemical 
potentials  of  relevant  species  in  the  membrane  are  bounded  by 
the  values  in  the  adjacent  reservoirs  (electrodes)  [12].  As  will  be 
demonstrated  here,  the  intuitive  notion  of  the  chemical  potentials 
of  species  in  the  membrane  being  bounded  by  the  values  in  the 
adjacent  reservoirs  (electrodes)  is  not  of  general  validity.  There  are 
conditions  under  which  chemical  potentials  in  the  membrane  may 
lie  outside  the  range  covered  by  the  two  electrodes  (reservoirs). 


Question:  How  does  Ms 
vary  through  the 
membrane? 


03 

c 

2 

_Q 

E 


x 


x  =  x0 


x  =  xf 


Enlarged  View  of  Interfaces 


Fig.  1.  A  schematic  of  a  membrane  separating  two  electrodes  (reservoirs).  The 
chemical  potentials  of  electrically  neutral  species  s,  /is,  in  the  two  electrodes  are 
p}s  and  /xf .  The  question  addressed  is:  how  does  /xs  vary  through  the  membrane? 
Positions  x0  and  Xf  are  in  the  reservoirs  (electrodes),  adjacent  to  the  membrane. 


Such  situations  can  lead  to  instability  of  the  membrane  and  thus  of 
the  device  in  question. 

The  following  discussion  begins  with  the  description  of  coupled 
processes  (or  chemical  reactions),  such  as  those  described  in  text¬ 
books  on  non-equilibrium  thermodynamics  (e.g.  Kondepudi  and 
Prigogine  [9],  Hasse  [10],  Denbigh  [11]).  Suppose  we  consider  two 
coupled  processes  (reactions).  One  of  them  occurs  such  that  its  free 
energy  decreases  as  the  process  occurs.  The  second  process  is  cou¬ 
pled  to  the  first.  However,  it  occurs  in  a  direction  such  that  its  free 
energy  increases  along  the  direction  of  reaction.  This  is  perfectly 
reasonable  as  long  as  the  combined  two-reaction  process  leads  to 
a  net  decrease  in  free  energy  [9].  What  does  this  mean  for  trans¬ 
port?  It  means  the  chemical  potentials  of  some  species  in  a  coupled 
transport  process  may  increase  along  the  direction  of  transport.  Our 
interest  is  in  determining  the  conditions  under  which  the  chemi¬ 
cal  potentials  of  one  or  more  electrically  neutral  species  exceed 
the  thermodynamic  stability  range  of  the  membrane.  We  now  con¬ 
sider  one  simple  but  an  illustrative  example  based  on  Stokes  law  of 
spheres  moving  though  a  fluid  to  demonstrate  the  role  of  coupling. 
Subsequently,  an  analogy  between  ion-electron  transport  through 
membranes  and  the  Stokes  law  problem  will  be  discussed.  The  illus¬ 
trative  example  of  two  spheres  immersed  in  a  fluid  (Stokes  law) 
describes  a  case  in  which  the  two  processes  are  fully  coupled. 


2.  Theory 

2.1.  Coupled  motion  of  spheres  in  a  fluid:  Stokes  law 

Fig.  2(a)  shows  two  spheres  of  radii  r\  and  r2  immersed  in  a 
fluid  of  viscosity  p,  joined  by  a  mass-less,  infinitely  thin  string  going 
over  a  frictionless  pulley  (pulley  outside  the  fluid).  Fig.  2(b)  shows 
two  spheres  immersed  in  a  viscous  fluid  connected  by  a  string,  but 
without  a  pulley.  As  given  in  Appendix  A,  the  velocities  of  the  two 
spheres  in  steady  state  (terminal  velocities  achieved)  with  a  string 
going  over  the  pulley  (Fig.  2(a))  are  given  by  [Eqs.  (Al)  and  (A2) 
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Pulley 


Fig.  2.  (a)  Two  spheres  of  masses  (corrected  for  buoyancy)  mi  and  m2,  of  radii  r\ 
and  r2  immersed  in  a  fluid,  connected  by  a  string  going  over  a  pulley  (pulley  outside 
the  fluid),  (b)  Two  spheres  of  masses  (corrected  for  buoyancy)  mi  and  m2,  of  radii 
ri  and  r2  immersed  in  a  fluid,  connected  by  a  string  but  without  a  pulley. 


/<,  molecule/atom/ion/electron,  at  some  position  r.  With  the 
passage  of  time,  ‘new’  k  molecules/atoms/ions/electrons  arrive 
at  r  and  force  Xk  =  -W  fik  refers  to  that  acting  on  new  k 
molecules/atoms/ions/electrons  at  r .  The  illustrative  example  given 
here  is  only  meant  to  demonstrate  that  the  coupling  terms  l\2X\ 
and  L2iX2  can  be  positive  or  negative,  and  this  determines  if 
the  chemical  potentials  of  some  species  increase  along  the  direc¬ 
tion  of  transport  during  coupled  transport.  We  will  now  examine 
atomic/ionic/electronic  transport. 

2.2.  Transport  through  predominantly  ionic  conductors 

Transport  of  electrically  charged  species  occurs  under  elec¬ 
trochemical  potential  gradients.  Even  in  predominantly  ionic 
conductors,  as  discussed  in  Appendix  B,  local  equilibrium  demands 
that  electronic  conductivity  while  may  be  negligible,  is  not  math¬ 
ematically  zero.  For  a  material  with  mobile  ions  i  of  valence  zit 
with  the  corresponding  neutral  species  being  s,  local  equilibrium 
at  r  between  the  neutral  species  s,  the  corresponding  constituent 
mobile  ion  i  of  valence  z,  ,  and  electrons,  e,  is  given  by 

s(r)  i(r)  +  Zj-e(r)  (24) 


from  Appendix  A] 

Vl  ~  V~6 3r*?(r1+r2)^  m 6nr](n  +  r2) 
=  -f'liV/Xi  -L12V/x2  =l\\X\  +  l\2X2 
and 

V2  ~  V  ~  67H7(ri  +  r2/  ™ ^  +  6707^1  +  r2) 
=  -T21V/X1  -L22V/X2  =T2iXi  +^22X2 


-m2g) 


-mg) 


(18) 


(19) 


The  coupling  coefficients  are  negative,  that  is 
1 


f  12  =  £21  =  -7 


<  0 


(20) 


67H7(ri  +r2) 

This  case  corresponds  to Ln  =L22  =  -Lu  =  -L2\  =  L=  l/(67H7(ri  +r2)) 
>  0.  For  the  case  without  a  pulley  (Fig.  2(b)),  the  velocities  are  [Eqs. 
(A7)  and  (A8)  from  Appendix  A] 

1  1 


V\  =  -v 


r(-mig)+  ; 


67H7(ri  +r2)  67777^1  +r2) 

=  -T11V/X1  -L12V/X2  =  TnXi  +L12X2 


(-mg) 


(21) 


and 

v2  =  -v  =  - — — ^ - s(-mg)  +  ^ — 7— - -S-mg) 

67H7(ri  +r2)  67H7(ri  +r2) 

=  -^2iV/Xi  -L22V/X2  —  L2\X]  +  L22X2  (22) 


Note  the  coupling  coefficients  are  positive,  that  is 

Li  2  —  L2I  =  — - 7 - r  >  0 

67H7(ri  +r2) 


(23) 


The  signs  of  the  forces  are  the  same  in  both  cases,  but  the  signs  of 
L12  =L2i  are  different.  Thus  in  the  two  cases  the  signs  of  the  terms 
l\2X2  and  L2\X\  are  different.  In  one  case,  in  which  L12X2  >  0  and 
L21X1  >0,  an  increase  in  energy  of  one  sphere  occurs  as  a  result  of 
coupling. 

The  example  of  spheres  joined  by  a  string  moving  in  a  fluid 
under  gravity  gives  a  simple  analogy  to  coupling.  In  most  sys¬ 
tems  of  interest,  transport  is  of  molecular/atomic/ionic/electronic 
species.  The  origin  of  coupling  in  transport  is  typically  Coulombic. 
Also,  processes  may  or  may  not  be  fully  coupled.  Addition¬ 
ally,  force  Xk  =  —  V  /xk  refers  to  that  acting  on  a  generic  species 


Analysis  of  transport  under  local  equilibrium,  the  fundamental 
assumption  in  virtually  all  transport  theory,  gives  an  equation  for 
the  chemical  potential  of  the  neutral  species  s  within  the  membrane 
in  terms  ionic  current  density  through  membrane,  /,-,  electronic  cur¬ 
rent  density  through  the  membrane,  Ie ,  and  ionic  and  electronic 
transport  resistances,  r((x)  and  r'e(x),  is  given  by  [Appendix  B,  Eq. 
(B32)] 


Ms(x)  =  MsO'o  )-zier'j(x)Ij+zier,e(x)le 
=  fJ'sCx'o)  ~ zie  |r-(x)/i  -  r'e(x)Ie J 


(25) 


which  is  in  a  form  similar  to  the  author’s  previous  derivation  [13]. 

Eq.  (25)  is  of  fundamental  interest  for  materials  exhibiting  pre¬ 
dominantly  ionic  transport  with  one  ionic  species  transporting, 
typical  of  electrolytes  with  only  one  ion  exhibiting  significant  con¬ 
ductivity. 

Returning  to  Fig.  1  letx0  refer  to  a  point  in  the  left  reservoir  (elec¬ 
trode),  close  to  the  membrane.  Let  Xf  refer  to  a  point  just  outside 
the  membrane  on  the  other  side  (in  the  right  reservoir/electrode 
close  to  the  membrane).  The  chemical  potentials  of  s  in  the  two 
reservoirs  (electrodes)  adjacent  to  the  membrane  are  pLs(x0)  and 
/xs(Xf).  These  are  the  experimentally  imposed  conditions;  for  exam¬ 
ple  in  a  lithium  ion  battery,  /xLi(x0)  and  /xLi(Xf)  could  be  the  lithium 
chemical  potentials  in  the  anode  and  the  cathode,  respectively.  The 
interface  thicknesses  are  assumed  to  be  very  small.  Thus,  properties 
and  thermodynamic  functions  vary  relatively  smoothly  through 
both  reservoirs  (electrodes)  and  the  membrane.  However,  they  vary 
sharply  across  the  ‘interfaces’  [  1 3,1 4].  As  properties  do  vary  through 
the  interfaces,  however  sharply,  one  can  still  define  effective  or 
average  interface  properties.  For  example,  area  specific  resistance 
to  ion  transport  across  the  interface  at  x0  may  be  given  by 


(26) 


where  or*0  is  the  average  ‘ionic  conductivity’  of  the  interface  at  x0 
and  a*,  is  the  thickness  of  the  interface  at  x0.  Eq.  (26)  also  gives  the 
interface  area  specific  resistance  as  an  integral,  if  such  a  variation 
through  the  interface  is  known.  This  information  on  spatial  varia¬ 
tion  through  the  interface  is  not  available,  and  the  integral  in  Eq. 
(26)  is  given  only  for  completeness.  The  parameters,  o*°  and  Sx0, 
are  also  not  separately  measurable.  However,  r?°  is  measurable. 
The  parameter  r*°  is  the  area  specific  charge  transfer  resistance 
and  thus  is  inversely  proportional  to  the  exchange  current  density 
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(excluding  any  effects  of  porous  electrode  contribution  to  electro¬ 
catalysis).  Similarly,  area  specific  resistance  to  ion  transport  across 
the  interface  at  Xf  is  rff .  Some  transport  of  electrons  must  also  occur 
across  both  interfaces.  The  corresponding  area  specific  resistances 
for  electron  transport  across  the  two  interfaces  are  given  similarly 
by  rx°  and  r*f. 

The  expectation  that  some  electron  transport  must  occur  across 
all  interfaces  (and  through  the  bulk)  can  be  understood  as  follows. 
All  materials  at  a  finite,  nonzero  temperature  exhibit  nonzero  free 
electron  concentration,  which  means  electronic  conductivity  can¬ 
not  be  mathematically  zero.  Thus,  some  transport  of  electrons  must 
occur  between  two  points  when  nonzero  electrochemical  potential 
gradient  of  electrons  is  established.  And  therefore  r*°  and  r*f  may  be 
very  large  compared  to  ion  transport  resistances,  but  are  not  math¬ 
ematically  infinite  for  reasons  of  local  equilibrium.  If  the  electronic 
resistances  were  infinite,  the  electronic  current  will  be  zero,  and 
from  Eq.  (25),  there  will  be  terms  of  the  type  ooxO  which  are  math¬ 
ematically  indeterminate,  which  means  the  chemical  potential  is 
undefined— violation  of  the  local  equilibrium  criterion  [13,15].  The 
very  low  level  electronic  conductivity  in  typical  ionic  conductors 
will  likely  be  of  hopping  or  a  polaron  type. 

The  chemical  potential  of  s  at  x0  (in  the  reservoir  or  the  elec¬ 
trode)  is  /zs(x0).  Consider  position  x  in  the  membrane  such  that 
x  >x0.  The  chemical  potential  of  s  at  x  is  given  by  Eq.  (25)  in  terms  of 
/zs(x0)  (with  x'0  replaced  by  x0)  in  which  the  net  area  specific  ionic 
and  electronic  resistances  are  respectively 


rf(x)  =  rf°  + 

and 

re(x)  =  r*°  + 


(X-Xo) 

(27) 

*1 

o 

* 

1  b 

(28) 

Eqs.  (27)  and  (28)  are  applicable  for  x0<x<Xf  and  the  /zs(x)  is 
given  by  Eq.  (25)  (with  x'0  replaced  by  x0).  In  order  to  relate  the 
chemical  potential  of  s,  /zs(xj),  at  Xj,  to  /zs(x0),  both  interfaces  must 
be  taken  into  account.  The  corresponding  area  specific  ionic  and 
electronic  resistances  are 


n(Xf)  =  r?°  + 

and 

re(xf )  =  rxe°  + 


(X/  —  x0) 

— - +  rff 

(29) 

O',  1 

(X/-X0) 

— - +  r^f 

Oe 

(30) 

The  corresponding  /rs(x/)  is  given  using  Eq.  (25 )  (with  x'0  replaced 
by  x0 )  by 


Ms(x/)  =  tis(Xo)  -  z,e  [rj(Xf)/,-  -  re(x/)JeJ  (31 ) 

The  /xs(x0)  and  /zs(xj)  are  the  experimentally  imposed  conditions 
(in  the  electrodes  or  the  reservoirs).  Our  interest  is  in  determining 
the  chemical  potential  of  s,  /zs(x)  within  the  membrane,  that  is  for 
x0<x<Xf  (Fig.  1).  Suppose  the  chemical  potentials  of  s  in  the  two 
reservoirs  (electrodes)  adjacent  to  the  membrane,  fis{x0)  and  /zs(xj), 
correspond  to  stable  conditions  for  the  membrane.  An  intuitively 
obvious  conclusion  that  therefore  the  membrane  must  be  stable 
‘under  transport’  is  not  of  general  validity,  as  discussed  in  what 
follows. 

Let  us  return  to  Eq.  (25)  (with  x'0  replaced  by  x0)  and  choose  the 
range  x0<x<Xf.  We  will  consider  two  cases. 


2.2.1.  Case  I 

The  ionic  and  the  electronic  currents  are  in  opposite  directions: 
That  is,  the  signs  of//  and  Ie  are  opposite.  The  term  [rz(x)/z  -  re(x)Je]  in 
Eq.  (25)  then  can  be  written  as  ±[r2-(x)  |/z|  +re(x)  |/e  |  ].  Thus  the /zs(x) 
for  x0  <  x  <  Xf  must  be  mathematically  bounded  by  /xs(x0 )  and  /xs(xj). 
That  is,  /zs(x0)  <  /xs(x)  <  /zs(xj)  or  /zs(x0)  >  /zs(x)  >  /zs(xj),  depending 


upon  whether  /zs(x0)  <  /zs(Xf)  or  /xs(x0)  >  /zs(Xf).  In  such  a  case,  the 
membrane  is  thermodynamically  stable  under  transport,  if  it  is  sta¬ 
ble  in  the  range  from  /zs(x0)  to  /zs(xj). 

2.2.2.  Case  II 

The  ionic  and  the  electronic  currents  are  in  the  same  direc¬ 
tion:  That  is,  the  signs  of  f  and  Ie  are  the  same  and  the  term 
[rz(x)/j  —  re(x)/e]  may  be  written  as  ±[q(x)  |/j|  -  re(x)  |/e|].  The  rel¬ 
ative  magnitudes  of  rz(x)  |/j|  and  re(x)  |/e|  now  determine  the  sign 
of  [rz(x)  |/j  |  -  re(x)  |  Je  |  ].  Thus  the  /zs(x)  need  not  be  mathematically 
bounded  by  /zs(x0)  and  /zs(xj).  Suppose  we  choose  /zs(x0)>/zs(Xf). 
Then,  the  preceding  implies  that  the  /zs(x)  for  x0<x<Xf  may  be 
higher  than  /zs(x0 ),  or  lower  than  /zs(Xf),  or  may  lie  in  the  range  cov¬ 
ered  by  fjLs(Xf)  and  /zs(x0).  If  /zs(x)  <  /zs(x/)  and/or  /zs(x)  >  /zs(x0),  and 
if  the  /xs(x)  lies  outside  the  thermodynamic  stability  range,  mem¬ 
brane  degradation  may  occur  even  when  exposed  to  reservoirs 
(electrodes)  corresponding  to  stable  conditions  for  the  membrane. 
Thus,  when  the  f  and  Ie  are  in  the  same  direction,  the  membrane 
may  become  thermodynamically  unstable  under  transport,  even 
when  exposed  to  stable  conditions  at  the  two  electrodes  (reser¬ 
voirs)  (/Xs(X0)  to  jL6s(X/)). 

2.3.  Abrupt  changes  across  interfaces 

It  is  commonly  assumed  in  much  of  the  electrochemical  litera¬ 
ture  that  reversible  electrodes  means  equilibration  of  the  chemical 
potential  of  the  neutral  species  (corresponding  to  the  mobile  ion  in 
the  electrolyte)  across  electrode/electrolyte  interfaces.  However, 
as  shown  in  Appendix  C,  this  is  generally  not  true.  The  transport 
theory  fundamentals  based  on  local  equilibrium  show  that  chem¬ 
ical  potentials  of  neutral  species  generally  exhibit  abrupt  changes 
across  interfaces,  even  when  electrodes  are  reversible.  In  general, 
abrupt  changes  in  chemical  potential  of  neutral  species  across  inter¬ 
faces  can  be  given  in  terms  of  ionic  and  electronic  current  densities, 
and  ionic  and  electronic  transport  parameters  of  the  interfaces 
[Appendix  C,  Eqs.  (C3)  and  (C4)]. 

2.4.  Implications  of  the  analysis  concerning  failure  of  lithium 
metal  and  li-ion  batteries 

Implications  of  the  above  analysis  concerning  lithium  metal  and 
lithium  ion  batteries  are  discussed  in  what  follows.  It  is  well  known 
that  lithium  batteries  (as  do  other  batteries)  degrade  under  various 
conditions.  Much  of  the  work  on  lithium  metal  and  lithium-ion  bat¬ 
teries  describes  roles  of  electrode  composition  and  microstructure, 
which  are  important  factors.  There  is  also  considerable  literature 
concerning  the  electrolyte,  especially  related  to  its  reactivity  with 
lithium  and  physical  properties  such  as  viscosity.  Most  electrolytes 
are  unstable  in  contact  with  a  lithium  metal  containing  anode  (and 
sometimes  may  also  be  unstable  in  contact  with  some  cathode 
materials)  and  form  a  solid  electrolyte  interphase  (SEI)  layer.  There 
is  extensive  literature  on  the  formation,  modeling  and  characteri¬ 
zation  of  the  SEI  layer  [  1 6-26],  and  its  role  in  dictating  cycle  life  and 
durability  [27].  The  present  work  suggests  that  transport  properties 
of  the  membrane,  including  the  SEI  layer  and  the  interfaces,  may 
be  important  in  cell  durability  in  light  of  the  model  presented  here. 

Lithium  metal  batteries  often  form  lithium  dendrites  during 
charging  which  can  short  the  battery.  Since  lithium  exists  in  a 
metallic  state  in  lithium  metal  battery  anodes,  the  likely  location  for 
dendrite  initiation  is  the  anode/electrolyte  interface.  This  is  exper¬ 
imentally  observed.  In  lithium  ion  batteries,  the  anode  is  typically 
graphite  in  which  lithium  is  intercalated  [28].  Thus,  the  chemical 
potential  of  lithium  in  the  anode,  is  lower  than  that  of  pure 
lithium,  and  this  should  lower  propensity  to  dendrite  forma¬ 
tion.  These  expectations  have  in  general  been  well  supported  by 
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observations.  Even  so,  lithium  metal  dendrites  or  lithium  plating 
are  known  to  occur  in  lithium  ion  batteries  under  certain  condi¬ 
tions.  These  observations  are  in  accord  with  the  analysis  presented 
here  which  shows  that  even  when  the  chemical  potential  of  lithium, 
/zLi,  in  both  electrodes  may  be  less  than  /zj^,  no  guarantee  exists 
that  it  would  so  remain  in  the  electrolyte  under  ah  conditions. 
The  model  presented  here  based  on  linear  non  equilibrium  ther¬ 
modynamics  shows  that  /zLi  in  the  electrolyte  may  exceed  that 
corresponding  to  pure  metallic  lithium,  /z°j,  causing  precipitation 
of  Li  as  dendrites  or  plating,  and/or  lead  to  electrolyte  destabiliza¬ 
tion  by  chemically  reacting  with  lithium.  This  scenario  is  possible 
if  under  any  conditions  lithium  ion  and  electronic  currents  through 
the  electrolyte  are  in  the  same  direction.  This  is  discussed  in  what 
follows. 

2.5.  Discharging  and  charging  ofli-ion  batteries 

Fig.  3(a)  is  a  schematic  of  a  Li-ion  cell  during  discharge.  Fig.  3(b)  is 
a  schematic  during  the  charging  step.  The  simple  schematics  do  not 
show  the  SEI  layers,  nor  do  they  show  the  electrode  details,  which 
are  assumed  to  be  of  uniform  composition  (ignores  concentration 
polarization  [29]).  The  objective  here  is  to  first  present  a  high-level 
description  of  the  differences  between  charge  and  discharge  steps 
insofar  as  the  chemical  potential  of  neutral  lithium,  /zLi(x),  at  x 
within  the  electrolyte  is  concerned.  The  x0  refers  to  the  anode,  close 
to  the  electrolyte.  The  Xf  refers  to  the  cathode,  close  to  the  elec¬ 
trolyte.  The  chemical  potential  of  lithium  in  the  anode  is  /z^  and 
that  in  the  cathode  is  /z^.  For  a  lithium  ion  battery,  /z^  <  /z^  <  /z£r 
With  reference  to  Fig.  3,  note  /z^  =  /zLi(x0)  and  /z^  =  /zLi(Xf ). 


During  discharge  (Fig.  3(a)),  lithium  ions  transport  from  the 
anode,  through  the  electrolyte,  to  the  cathode.  Thus,  the  current 
due  to  lithium  ion  transport,  /Li  =  f  is  from  the  anode  through  the 
electrolyte  to  the  cathode.  The  cathode  is  at  a  higher  electric  poten¬ 
tial,  (p}  than  the  anode,  that  is  cpc  >  (fA  and  electron  transport  through 
the  electrolyte ,  however  small  may  be  the  electronic  conductivity  in 
relation  to  the  ionic  conductivity,  occurs  from  the  anode  through 
the  electrolyte  to  the  cathode.  Thus,  the  electronic  current  through 
the  electrolyte,  Je,  is  from  the  cathode  through  the  electrolyte  to  the 
anode  since  electrons  are  negatively  charged.  Therefore,  during  dis¬ 
charge,  h  and  Ie  are  in  opposite  directions,  and  according  Eq.  (25), 
/zLi(x)  in  the  electrolyte  is  bounded  by  the  values  at  the  electrodes, 
that  is  ficu  <  /zLi(x)  <  /z^.  In  such  a  case,  lithium  precipitation  in 
the  electrolyte  cannot  occur  since  /z^  <  jz°  . 

During  charging  as  shown  in  Fig.  3(b),  lithium  ions  transport 
from  the  cathode,  through  the  electrolyte,  to  the  anode.  Thus,  the 
ionic  current  f  is  from  the  cathode  through  the  electrolyte  to  the 
anode.  The  cathode  is  still  at  a  higher  electrical  potential  than  the 
anode  (even  more  so),  that  is  (pc  >  ((A  and  electronic  current  through 
the  electrolyte  Ie  is  still  from  the  cathode  to  the  anode.  Thus,  during 
charging,  f  and  Ie  are  in  the  same  direction.  The  relative  magnitudes 
of  r,-(x)/j  and  re(x)/e  now  determine  the  sign  of  [r^x)/,-  -re(x)/e]  in 
Eq.  (25)  (with  x'0  replaced  by  x0)  and  determine  /zLi(x).  Thus  dur¬ 
ing  the  charging  step,  /zLi(x)  in  the  electrolyte  need  no  longer  be 
bounded  by  the  values  at  the  electrodes  (by  /z^  and  /z^ ).  Conditions 
could  develop  under  which  the  /zLi(x)  may  exceed  /z^  and  this  will 
lead  to  the  precipitation  of  lithium  metal  in  the  electrolyte.  This 
precipitation  would  likely  occur  at  the  anode/electrolyte  interface 
and  would  appear  as  a  dendrite  or  as  plating.  An  increase  in  /zLi(x) 
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Fig.  3.  (a)  A  schematic  of  a  lithium-ion  cell  during  discharge.  The  directions  of  the  ionic  current,  /*,  and  the  electronic  current,  Ie,  through  the  electrolyte  are  shown.  The  ionic 
and  the  electronic  currents  through  the  electrolyte  are  in  opposite  directions.  Note  that  |L|«|  /,  | .  (b)  A  schematic  of  a  lithium-ion  cell  during  charge.  The  directions  of  the 
ionic  current,  /,-,  and  the  electronic  current,  Je,  through  the  electrolyte  are  shown.  The  ionic  and  the  electronic  currents  through  the  electrolyte  are  in  the  same  direction.  Note 
that  \Ie\  «;  |/,-|. 
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/i  =  A  +  4=|/,|+|/e|=/,  v|/r|«|/,| 

Fig.  4.  A  schematic  of  a  lithium-ion  cell  with  an  SEI  layer  at  the  anode  during  charge. 
Note  that  \Ie\  «;  |  /,■  | . 


Ionic  and  electronic  conductivities  of  the  SEI  layer  are  respectively 
trp1  and  ofEI,  and  its  thickness  is  <$sei-  The  equivalent  circuit  thus 
contains  four  series  segments:  (1 )  the  anode/SEI  layer  interface,  (2) 
the  SEI  layer,  (3)  the  electrolyte,  and  (4)  the  electrolyte/cathode 
interface.  Each  segment  has  two  legs  in  parallel;  an  ionic  leg  and 
an  electronic  leg.  Each  ionic  leg  contains  an  internal  EMF,  which  is 
a  measure  of  the  difference  in  /xy  across  that  segment.  Thus,  the 
internal  EMFs  are 
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F*„  <  ~  lh° 

1  e 

(32) 

Sei  Mu  A^Li 

Ei  =  e 

(33) 

£el 

Ei  =  e 

(34) 

and 

m  ^  -  '4 

1  e 

(35) 

during  transport  may  also  lead  to  a  reaction  of  lithium  with  the 
electrolyte  causing  electrolyte  degradation.  The  reason  anode  SEI 
layer  forms  is  because  the  electrolyte  is  unstable  at  high  thermo¬ 
dynamic  lithium  activity  (chemical  potential).  If  /xLi(x)  increases  to 
a  value  at  which  reaction  can  occur,  it  most  likely  will  occur  since 
the  electrolyte  is  typically  a  liquid  or  a  gel,  and  the  kinetics  of  such 
reactions  will  be  fast.  This  will  lead  to  internal  destabilization  of 
the  electrolyte  and/or  further  thickening  of  the  SEI  layer. 

2.6.  The  role  of  SEI  layer  in  preventing/suppressing  degradation 

An  SEI  layer  may  form  at  both  the  anode  and  the  cathode.  Fig.  4 
shows  a  schematic  of  a  cell  with  an  SEI  layer  at  the  anode,  which 
may  not  be  of  uniform  composition,  nor  of  uniform  morphology. 
The  SEI  layer  here,  however,  is  assumed  to  be  of  uniform  prop¬ 
erties.  Also,  it  is  assumed  that  the  SEI  layer  forms  over  the  entire 
anode/electrolyte  interface,  is  of  uniform  thickness  and  anode  exfo¬ 
liation  does  not  occur.  The  SEI  layer  is  predominantly  an  ionic 
(lithium  ion)  conductor.  Its  lithium  ion  conductivity  albeit  is  lower 
than  of  the  electrolyte.  But  we  expect  its  electronic  conductivity  to 
be  very  low  [16].  The  rate  at  which  an  SEI  layer  thickens  over  time 
(even  when  the  cell  is  not  in  use)  is  a  measure  of  lithium  permeation 
from  the  anode  through  the  SEI  layer  [30].  This  means  electronic 
conduction  through  the  SEI  layer  may  be  small,  but  is  not  mathe¬ 
matically  zero.  In  what  follows,  we  will  examine  conditions  under 
which  the  /xLi(x)  in  the  electrolyte  equals  or  exceeds  /x° .  Since  this 
may  occur  during  charging,  this  will  be  the  case  examined  first. 
A  convenient  approach  is  to  describe  the  cell  using  a  simplified 
equivalent  circuit,  as  described  below. 


X' 

/xL?  is  the  /xLi  in  the  SEI  layer  very  close  to  the  anode/SEI  layer 
interface  and  /x£j  is  the  /xLi  in  the  electrolyte  very  close  to  the 
electrolyte/cathode  interface.  The  Nernst  EMF  of  the  cell  is  given 
by 


A  C 

En  =  Ef»  +  Efa  +  EfL  +  Eff  =  Mu  ~  (36) 

The  ionic  and  the  electronic  area  specific  resistances  of  the  SEI 
layer  are 


<*>SEI 

CTp1 


and  rp1 


^SEI 

ofH 


(37) 


The  total  area  specific  ionic  and  electronic  resistances  of  the  cell 
are  respectively 


R:  =  rf°  +  rp1  +  rp  +  rff 


(38) 


and 


Re  =  r*°  +  rp1  +  rp  +  r*f  (39) 

Consider  now  charging  of  the  cell.  The  application  of  Eq.  (25) 
(with  x'0  replaced  by  x0)  for  /xLi  at  a  distance  x-x0  where  x  is  a 
position  in  the  electrolyte  gives 


Mii(x)  =  Mu  +  e 


^o+rSE.+  (*  —  |,.| 


(x  -x0)  -  (SsElA 
cre  J 


Ie 


(40) 


2.7.  Equivalent  circuit 

Assuming  transport  properties  are  uniform  within  a  given  layer, 
and  assuming  no  interfacial  losses  occur  at  the  SEI  layer/electrolyte 
interface,  a  Li-ion  cell  operating  in  steady  state  may  be  represented 
by  a  simplified  equivalent  circuit  shown  in  Fig.  5.  The  equivalent  cir¬ 
cuit  given  here  only  represents  the  region  from  x0  to  Xj,  and  ignores 
both  electrodes  but  includes  the  interfaces.  Both  electrodes  are 
assumed  to  be  of  uniform  (but  different)  compositions  and  at  con¬ 
stant  (but  different)  electric  potentials,  (fA  and  qf.  If  they  are  not  of 
uniform  composition  and  potentials  (when  concentration  polariza¬ 
tion  effects  are  significant  [29]),  then  /x^  and  qA  refer  to  the  values 
in  the  anode  close  to  the  anode/SEI  interface,  and  /x^  and  qf  refer 
to  the  values  in  the  cathode  close  to  cathode/electrolyte  interface. 


Let  us  set  x  -  x0  =  <$Sei,  which  corresponds  to  the  SEI/EL  interface. 
The  chemical  potential  of  lithium  at  the  SEI/EL  interface  is  then 
given  by 

Mu(SEI/EL)  =  +  e  [(r?°  +  r?EI )  |/,|  -  (tf>  +reSEI)  |/e|J  (41) 

Now, 


Ea  ~  En 

Ri 


and 


Ea  -  En 


Ea  _ _ Ea _ 

Re  r*°  +  rp1  +  rp  +  r*f 


(42) 


(43) 
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Fig.  5.  An  equivalent  circuit  during  charging  in  steady  state  for  the  cell  schematic  shown  in  Fig.  4.  The  equivalent  circuit  contains  four  segments:  (a)  anode/SEI  layer  interface, 
(b)  SEI  layer,  (c)  electrolyte,  and  (d)  electrolyte/cathode  interface.  Each  segment  has  an  ionic  leg  and  an  electronic  leg  in  parallel.  The  circuit  elements  in  between  the  filled 
circles  are  not  physically  separable.  Note  that  \Ie  |  «  \h\. 


Thus,  substituting  in  Eq.  (40),  the  chemical  potential  of  lithium 
at  the  SEI/EL  interface,  /xLi(SEI/EL),  is  given  by 


/xLi(SEI/EL)  =  ,iAu+e 


(r?°  +  r?EI)  (Ea  -  En ) 

^j~Xo  _|_  j-SEI  _|_  j.EL  _|_  pXj  ^ 


(r?  +  r|EI)  Ea 


(r£°  + 


.SEI  I  rEL 
e  '  1  e 


+  r*ef) 


(44) 


Eq.  (44)  gives  the  chemical  potential  of  lithium  at  the  SEI/EL 
interface  during  charging  in  terms  of  the  Nernst  voltage  {EN),  the 
applied  voltage  ( EA ),  the  chemical  potential  of  lithium  in  the  anode 
(/x^)  and  the  various  transport  resistances. 


2.8.  Four  cases 


Many  different  scenarios  are  possible.  We  will  now  examine  four 
cases.  These  are  selected  to  illustrate  the  application  of  Eq.  (44)  and 
the  resulting  implications. 


2.8.1.  Casel 

Suppose  the  SEI  layer  exhibits  very  high  ionic  and  electronic 
resistances  such  that 


r.  »  rf°,  rEL,  rff  and  r; 


SEI 


»  r*°,  rfL,  r*s 


Then  Eq.  (44)  reduces  to 


Atu(SEI/EL)  «  Mu  -  eEN  «  Mu  (45) 

That  is  /iiiiSEIlEL )  is  about  the  same  as  in  the  cathode.  The  /xLl- 
at  the  anode/SEI  layer  interface  (in  the  SEI  layer)  is  also  the  same, 
that  is 

I*U  ~  Mu  -  eEN  «  Mu  (46) 

Fig.  6(a)  shows  the  corresponding  schematic  variation  of  /xLi(x) 
through  the  electrolyte  and  the  SEI  layer.  The  /xLi  in  the  electrolyte 
and  in  the  SEI  layer  are  about  the  same  as  in  the  cathode,  and  no  pre¬ 
cipitation  or  plating  of  metallic  lithium  is  expected  nor  is  a  reaction 
with  the  electrolyte  expected  (assuming  the  electrolyte  is  stable 
corresponding  to  the  cathode  chemical  potential  of  Li,  /x^).  This 


however  comes  with  a  penalty— the  cell  ionic  resistance  is  rela¬ 
tively  high  due  to  high  rSEI,  thus  lowering  the  performance  (current 
density  and  power  density). 

2.8.2.  Case  2 

Suppose  the  SEI  layer  is  a  good  ionic  conductor  but  a  poor  elec¬ 
tron  conductor.  Then  rfEI  »  r£°,  rEL,  r£f  and  Eq.  (44)  then  reduces 
to 


Mu(SEI/EL)  «  Mu  -  e 


(rf°  +  rSEI)  En  +  (rEL  +  rff )  EA 
(rf°  +  rSEI  +  rEL  +  rffj 


(47) 


The  term  in  square  brackets  of  Eq.  (47)  is  positive.  Thus, 
Mu(SEI/EL)  <  /x^  and  precipitation/plating  of  lithium  or  reaction 
with  the  electrolyte  is  not  expected.  The  preceding  implies  that  the 
SEI  layer  should  be  an  excellent  electronic  insulator.  This  of  course 
is  well  known.  However,  the  previously  discussed  reason  for  the 
desired  high  electronic  resistance  of  the  SEI  layer  is  to  prevent  the 
SEI  layer  from  growing  too  fast  which  can  shorten  cell  life  and  lower 
capacity.  The  present  analysis  shows  that  it  is  also  important  to 
lower  the  /xLi  in  the  electrolyte  and  lower  the  propensity  for  cell 
degradation  during  charging,  which  is  achieved  if  rfEI  is  large.  The 
present  work  also  suggests  that  anode  purity  may  be  an  important 
consideration.  Especially,  it  should  be  free  of  transition  metal  ions 
which  may  increase  the  electronic  conductivity  of  the  SEI  layer. 


2.8.3.  Case  3 

Now  suppose  rSEI  »  rXo ,  rEL,  rx/  but  now  rfEI,  rx°  rEL,  rx/.  This 
could  be  the  case  if  some  transition  metal  ions  are  present  in  the 
SEI  layer  or  some  metal  or  graphite  particles  inadvertently  pene¬ 
trate  the  SEI  layer  (due  to  defects  formed  during  manufacturing, 
such  as  a  rough  anode/electrolyte  interface),  or  impurities  in  the 
anode  which  could  introduce  localized  or  non-localized  electronic 
conductivity  in  the  SEI  layer.  Then  Eq.  (44)  reduces  to 

Mu(SEI/EL)  S3  Mu  +  e(EA  -  En)  =  Mu  +  ?EA  (48) 

During  charging,  we  must  have  EA>E^.  This  means 
Mu(SEI/EL)  >  Mu  (49) 

That  is,  now  /xy  in  the  electrolyte  is  not  bounded  by  the  val¬ 
ues  at  the  electrodes.  Typically,  lithium  potential  as  intercalated  in 
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Case  1  Case  3 


lL  =  l,  +  le=\l,\+\le\  =  li  v|4|«|/,| 

(a) 


/i  =  //+4=|/i|+|4|=/i 

(b) 


Fig.  6.  (a)  A  schematic  of  a  cell  for  Case  1 :  the  chemical  potential  of  lithium  in  the  electrolyte,  /xu(x),  is  about  the  same  as  at  the  cathode,  /x[..  In  such  a  case,  precipitation  of 
lithium  within  the  cell  is  not  expected.  Note  that  \Ie\  «;  |/,- 1 .  (b)  A  schematic  of  a  cell  for  Case  3:  the  chemical  potential  of  lithium  in  the  electrolyte  at  the  SEI/EL  interface 
can  exceed  that  at  the  anode,  and  may  even  reach/exceed  that  corresponding  to  pure  lithium,  In  such  a  case,  precipitation  of  lithium  within  the  cell  may  occur  leading 
to  degradation  of  the  cell.  Note  that  \Ie\  «;  \h\. 


graphite  is  ~0.25  V  with  respect  to  pure  lithium.  This  means 

—  ~  Ml'~°-25V  (5°) 

Thus,  if  the  applied  voltage  EA  exceeds  EN  by  0.25  V, 

Mu(SEI/EL)  >  /Zy  (51) 

If  this  situation  occurs,  lithium  metal  can  precipitate  within 
the  electrolyte,  initiating  at  the  SEI/EL  interface.  The  cell  may  then 
degrade  by  various  mechanisms— formation  of  lithium  dendrites 
or  plating,  electrolyte  degradation  (reaction  with  lithium  due  to 
high  chemical  potential  of  lithium)  leading  to  further  thickening  of 
the  SEI  layer  or  both.  Fig.  6(b)  shows  a  schematic  variation  of  /xLj(x) 
through  the  electrolyte  for  this  case.  The  preceding  thus  shows  that 
even  a  modest  level  of  electronic  conduction  in  the  SEI  layer  is  detri¬ 
mental  to  cell  stability  and  can  lead  to  dendrite  formation/plating 
during  charging. 

2.9.  The  proposed  mechanism  and  implications  concerning 
capacity  loss 

As  discussed  here,  during  discharge  (single  cell  only),  /xLi  in  the 
membrane  (electrolyte)  is  bounded  by  the  values  at  the  electrodes. 
Thus  high  values  of  /xLi  in  the  electrolyte  are  not  expected.  How¬ 
ever,  during  charging,  /xLi  in  the  electrolyte  may  exceed  /x^.  The 
most  likely  location  of  high  /xLi  will  be  at  (or  close  to)  the  SEI/EL 
interface.  Should  this  occur,  dendrites  (or  plating)  could  occur  or 
electrolyte  may  react  with  lithium  to  form  (further  thicken)  the 
SEI  layer,  or  both.  The  source  of  this  lithium  is  the  cathode  during 
charging.  During  successive  discharge-charge  cycles  less  and  less 
active  lithium  will  be  available,  and  the  step  during  which  capacity 
loss  occurs  will  primarily  be  the  charging  step. 

Recently,  Zhou  and  Notten  [31  ]  carried  out  an  interesting  study 
on  lithium-ion  cells.  They  subjected  a  number  of  lithium-ion  cells  to 
discharge-charge  cycles.  The  initial  capacity  was  ~1 00  mAh,  which 
degraded  to  ~80  mAh  at  the  end  of  the  test  protocol.  The  last  step 
for  each  cell  was  a  full  discharge.  They  then  separated  the  cell  com¬ 
ponents  and  new  cells  of  two  types  were  constructed.  Type  I  cell : 
These  cells  were  made  with  a  new  separator,  a  new  anode,  but  a 
cathode  from  the  cycled  cells.  Type  2  cell :  These  cells  were  made 


with  a  new  separator,  a  new  cathode,  but  an  anode  from  the  cycled 
cells.  Type  1  (old  cathode)  cells  had  a  capacity  of  ~80mAh,  the 
same  as  the  original  cells  at  the  end  of  cycling.  Type  2  (old  anode) 
cells,  by  contrast,  had  a  capacity  of  ~100mAh,  the  same  as  the 
pristine  cells.  These  results  thus  demonstrated  that  the  capacity 
loss  could  be  attributed  to  the  charging  step  during  which  loss  of 
lithium  occurred.  These  observations  are  in  accord  with  the  theory 
developed  here  showing  the  anticipated  increase  in  the  chemi¬ 
cal  potential  of  lithium  within  the  electrolyte  during  the  charging 
step  and  its  plating  or  reaction  with  the  electrolyte  causing  loss  of 
capacity,  the  origin  of  which  lies  in  the  Onsager  theory. 

2.9 A.  Case  4 

An  electronically  conductive  separator:  A  typical  separator  for 
a  lithium-ion  battery  is  a  micro-porous  polymeric  material  such  as 
polyethylene  or  polypropylene.  Both  of  these  materials  are  elec¬ 
tronic  insulators  [32].  A  lithium  salt  solution  in  a  non  aqueous, 
organic  polar  liquid  such  as  ethylene  carbonate  is  impregnated  into 
such  a  micro-porous  separator.  We  will  now  examine  the  possible 
use  of  a  micro-porous  separator  which  has  a  rather  high  electronic 
conductivity.  This  is  a  radically  different  approach  from  the  current 
practice  [32  ].  As  discussed  in  what  follows,  a  significant  benefit  may 
be  achieved  in  suppressing  degradation  and  capacity  fade  provided 
the  SEI  layer  is  an  excellent  electronic  insulator  but  the  separator 
has  a  relatively  high  electronic  conductivity. 

Consider  now  that  rEL  &  0  (the  separator  is  a  good  electronic 
conductor )  and  also  assume  r*f  «  0.  Eq.  (44)  then  becomes 

Mu(SEI/EL)  ~  My  +  e 
or 

/iLi(SEI/EL)  «=  /4  -  e 

The  term  in  the  square  brackets  in  Eq.  (53)  is  positive.  Thus, 
/xLi(SEI/EL)  <  /x^.  Note  Eq.  (53)  is  the  same  as  Eq.  (47).  The  dif¬ 
ference  is  that  Eq.  (47)  was  derived  assuming  the  SEI  layer  is  an 
electronic  insulator.  Eq.  (53)  is  derived  assuming  the  separator  is  a 
good  electronic  conductor.  Effectively,  both  achieve  the  same  result 


(rf°  +  rSEI)  En  +  (rEL  +  rx/)  EA 
(rf°  rSEI  _|_  rEL  _|_  rxf^j 


(53) 


(rf°  +  rp1)  (Ea  ~  En) 

_|_  ySEI  _|_  j-EL  -|-  ^ 
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Fig.  7.  A  schematic  of  a  lithium-ion  cell  with  a  bi-layer  separator. 


insofar  as  the  /zLi(SEI/EL)  is  concerned.  It  is  understood  of  course 
that  the  SEI  layer  must  exhibit  high  electronic  resistance  to  min¬ 
imize  internal  leakage.  What  Case  4  shows  is  that  the  lowering 
of  /zLi(x)  may  also  be  achieved  by  using  a  separator  of  high  elec¬ 
tronic  conductivity.  In  such  a  case,  precipitation  of  lithium  should 
not  occur  and  also  no  further  reaction  with  electrolyte  should 
occur  (except  with  whatever  lithium  that  permeates  from  the 
anode  through  the  SEI  layer).  Also,  then  the  SEI  layer  functions 
as  the  electrolyte  and  the  micro-porous  layer  with  the  electrolyte 
impregnated  in  it  is  a  composite  mixed  ionic  electronic  conduc¬ 
tor.  Now,  however,  because  of  high  electronic  conduction  through 
the  separator,  /zLi(x)  <  /z^  which  suppresses  dendrite  formation 
(or  plating).  The  present  analysis  thus  shows  that  the  separator 
actually  could  be  made  of  a  reasonably  good  electron  conducting 
material— and  need  not  be  an  electronic  insulator  as  has  been  the 
common  practice.  The  analysis  also  suggests  that  the  electrolyte 
and  the  anode  should  be  essentially  free  of  any  transition  met¬ 
als  to  minimize  electronic  conduction  through  the  SEI  layer.  An 
SEI  layer  free  of  transition  metals  should  be  a  very  good  electronic 
insulator  and,  with  high  electronic  conduction  occurring  through 
the  separator  will  maintain  /zLi(SEI/EL)  <  /z^  and  thereby  suppress 
degradation. 

2. 1 0.  Bi-layer  separator 

The  preceding  discussion  suggests  that  an  ideal  separator  may 
be  a  bi-layer  porous  structure  consisting  of  a  relatively  thick  porous 
electronic  conductor  and  a  thin  porous  electronic  insulator.  This 
could  be  achieved,  for  example,  by  applying  a  coating  of  polyethy¬ 
lene  or  polypropylene  on  a  porous  carbon  paper.  The  anode  should 
be  in  contact  with  the  electronically  insulating  layer  (polyethylene 
or  polypropylene)  while  the  cathode  should  be  in  contact  with  the 
electronically  conducting  layer  (carbon  paper).  In  this  manner,  the 
chances  of  the  electronic  conducting  layer  coming  in  direct  con¬ 
tact  with  the  anode,  in  the  case  of  breach,  will  be  minimized.  Fig.  7 
shows  a  schematic. 

2. 11.  Battery  pack  failure  due  to  cell  imbalance 

We  have  seen  that  for  a  single  cell  during  discharge,  /zLi(x)  is 
always  bounded  by  the  values  at  the  electrodes.  As  discussed  below, 
such  however  is  not  the  case  with  battery  packs  containing  many 
series-connected  cells.  The  following  discussion  is  restricted  to  dis¬ 
charging  of  a  lithium  ion  battery  containing  a  number  of  lithium-ion 
cells  connected  in  series.  In  a  typical  laptop  lithium-ion  battery, 


4-6  cells  are  connected  in  series.  For  higher  voltage  needs,  several 
cells  are  connected  in  series,  such  as  for  automotive  applications. 
In  the  Tesla  Roadstar  car  battery  the  nominal  operating  voltage  is 
~375  V,  indicating  over  1 00  cells  connected  in  series.  The  following 
discussion  is  particularly  relevant  to  such  applications. 

Consider  N  lithium-ion  cells  connected  in  series.  Let  us  assume 
that  N-  1  of  them  have  an  identical  internal  resistance,  R ,  but  one 
cell  has  a  much  higher  resistance,  R';  that  is  R'>R.  Suppose  the 
battery  is  discharged  at  a  voltage  V,  with  the  average  cell  voltage 
Eceii =  WN.  Since  the  cells  are  connected  in  series,  the  same  net  cur¬ 
rent,  /,  flows  through  each  cell.  Each  cell  also  has  the  same  Nernst 
voltage,  En.  The  battery  voltage  is  V=  NEcell  =  NEn  -  I(N  -  1  )R  -  IR'. 
The  voltage  across  each  of  the  N-  1  identical  cells  is  v  =  EN-IR 
while  the  voltage  across  the  high  resistance  cell  is  v'  =  EN-IR'. 
The  average  cell  voltage  is  Ecell  =EN  -J(((N-  \)IN)R  +  {R'  IN)).  For 
1,  the  voltage  across  each  of  the  N-  1  identical  cells  is  about 
v  «  Ece ii,  which  is  the  average  cell  voltage.  However,  that  across  the 
unbalanced  cell  is  v'  =  EN-  IR'.  If  R'  is  large,  the  voltage  across  this 
cell  may  become  negative,  that  is  v'  =  EN-  IR'  <  0.  This  scenario  is 
depicted  in  Fig.  8.  The  direction  of  the  ionic  current,  f  (the  mag¬ 
nitude  of  which  is  about  I)  is  the  same  through  all  cells.  However, 
the  unbalanced  cell  operating  at  a  negative  voltage  has  its  elec¬ 
tronic  current,  I'e ,  reversed  in  direction  compared  to  the  normal 
cells.  That  is,  in  the  unbalanced  cell,  the  directions  of  the  ionic  and 
the  electronic  currents  are  the  same.  In  this  cell,  the  /zLi(x),  need  no 
longer  be  bounded  by  the  values  at  the  electrodes  (by  /z^  and  /z^) 
and  lithium  metal  precipitation  may  occur  in  the  electrolyte.  For 
the  unbalanced  cell,  the  chemical  potential  of  lithium  at  x  in  the 
electrolyte  is 

MliM  =  Mu  -  e[/r;(x)  -  I'er'e(x)]  (54) 

where  r'(x)  is  the  ionic  area  specific  resistance  between  the  anode 
and  position  x  in  the  electrolyte,  r'(x)  is  the  electronic  area  spe¬ 
cific  resistance  between  the  anode  and  position  x  in  the  electrolyte, 
and  I'e  is  the  electronic  current  flowing  through  the  cell.  In  this 
unbalanced  cell,  the  anode  is  at  a  higher  electric  potential  than  the 
cathode.  That  is,  for  any  value  of  x  in  the  electrolyte,  qA’  >  <p(x)  > 
(pc .  Eq.  (54)  may  also  be  given  by 

MuW  =  Mu  +  e  [(</'  -  <PW)  -  /r!(x)J  (55) 

If  the  term  in  the  square  brackets  in  Eq.  (55)  is  positive,  the  /zLi(x) 
in  the  electrolyte  will  exceed  /z^.  If  the  /zLi(x)  exceeds  /z^,  lithium 
may  precipitate  in  the  electrolyte  or  react  with  it.  A  schematic  plot 
of  /zLi(x)  vs.  x  in  such  an  unbalanced  cell  is  shown  in  Fig.  9.  The 
expectation  is  that  high  /zLi  will  occur  in  the  electrolyte  close  to 
the  electrolyte/cathode  interface.  Thus,  in  such  a  case,  degrada¬ 
tion  may  initiate  at  the  cathode/electrolyte  interface.  Also,  an  SEI 
layer  due  to  high  /zLi  will  form  at  the  cathode  (Fig.  9)  in  addition 
to  that  at  the  anode.  The  present  work  thus  shows  that  there  are 
some  post  mortem  indicators  one  should  look  for  in  determining 
the  possible  reasons  for  battery  degradation/failure.  The  present 
work  suggests  that  in  a  battery  pack  unbalanced  cells  may  exhibit 
SEI  layers  at  both  anode/electrolyte  and  cathode/electrolyte  inter¬ 
faces.  Cells  behaving  normally  by  contrast  will  exhibit  SEI  layers 
only  at  anode/electrolyte  interfaces.  These  statements  assume  that 
the  electrolyte  under  normal  conditions  is  stable  with  respect  to  the 
cathode.  If,  however,  an  SEI  layer  also  forms  at  the  cathode  under 
normal  operating  conditions,  then  it  will  mean  that  an  unbalanced 
cell  in  a  battery  pack  will  exhibit  a  thicker  SEI  layer  at  the  cathode 
compared  to  the  normal  cells. 

2. 12.  A  comparison  with  the  Stokes  law  problem 

The  analysis  of  two  spheres  connected  by  a  string  immersed 
in  a  fluid  moving  under  gravity  was  conducted  for  two  cases 
(Appendix  A):  (a)  string  over  a  pulley,  and  (b)  no  pulley.  In  the 
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Fig.  8.  (a)  A  schematic  of  a  lithium-ion  battery  containing  N  cells  in  series  of  which  N-  1  are  identical  and  one  (the  unbalanced  cell)  has  a  resistance  much  larger  than  the 
other  cells,  (b)  A  schematic  of  voltage  vs.  current  plots  for  the  N  -  1  identical  cells  and  for  the  unbalanced  cell.  At  a  high  enough  current,  the  unbalanced  cell  may  switch 
polarities.  In  such  a  case,  the  ionic  and  the  electronic  currents  will  be  in  the  same  direction  in  the  unbalanced  cell  during  discharge  leading  to  cell  failure,  which  may  eventually 
lead  to  battery  failure. 


former,  the  lighter  sphere  moves  up  (moves  along  the  direction  of 
increasing  ‘chemical  potential’).  In  the  latter,  both  spheres  move 
down— decrease  their  ‘chemical  potentials’. 

The  transport  of  ionic,  electronic,  and  neutral  species  examined 
here  (Appendix  B)  is  analogous  to  these  two  cases.  Note  that  the 
transport  of  each  charged  species  was  assumed  to  occur  down  its 
electrochemical  potential  gradient.  This  is  a  common  assumption 
made  in  virtually  all  reported  studies.  This  assumption  naturally 
means  that  in  all  cases  considered  here,  ionic  and  electronic  species 
always  move  down  their  electrochemical  potential  gradients.  How- 


Unbalanced  Cell  in  a  Battery  Pack 
Discharge  Step 

SEI  Formation  at  Cathode 


N-1  Cells  +  Load 


Current 


Fig.  9.  A  schematic  showing  the  variation  of  lithium  chemical  potential,  /xLi(x), 
through  the  electrolyte  during  discharge  in  an  unbalanced  cell  in  a  battery.  The  /xLi(x) 
may  increase  with  x,  and  may  reach  a  high  value  in  the  electrolyte,  close  to  the  elec¬ 
trolyte/cathode  interface  If  exceeds  /x°,  lithium  may  precipitate  or  may 
react  with  the  electrolyte.  An  SEI  layer  then  may  also  form  at  the  cathode/electrolyte 
interface.  Note  that  \Ie\  «:  |/j  | . 


ever,  such  is  not  the  case  with  the  corresponding  electrically  neutral 
species— and  therein  lies  the  analogy. 

Consider  for  example  discharging  of  a  lithium-ion  cell  shown 
in  Fig.  3(a).  The  spatial  variation  of  /xLi(x)  as  a  function  position  x 
through  the  electrolyte  is  such  that  it  decreases  from  the  anode  to 
the  cathode.  This  case  is  analogous  to  spheres  joined  by  a  string, 
falling  under  gravity  (no  pulley).  All  species  move  down  their 
respective  thermodynamic  (chemical  or  electrochemical)  potential 
gradients. 

Now  consider  charging  of  a  lithium-ion  cell  shown  in  Fig.  6(b). 
Again,  lithium  ions  and  electrons  move  down  their  respective 
electrochemical  potential  gradients  (not  surprising— this  was  the 
assumption  made  in  the  problem  formulation,  a  standard  assump¬ 
tion  in  virtually  all  reported  studies  on  charged  species  transport). 
However,  note  that  electrically  neutral  lithium  is  effectively  trans¬ 
ported  up  its  chemical  potential  gradient  through  the  electrolyte. 
Similar  is  the  case  shown  in  Fig.  9  for  discharge  through  an  unbal¬ 
anced  cell  in  a  battery  pack.  These  cases  are  analogous  to  spheres 
connected  by  a  string  going  over  a  pulley— the  lighter  sphere  moves 
up  its  ‘chemical  potential’  gradient. 

3.  Summary 

A  theory  of  degradation  of  active  electrochemical  devices  based 
on  linear  non-equilibrium  thermodynamics  is  presented.  It  is 
shown  that  even  in  predominantly  ionic  conductors  electronic  con¬ 
duction  cannot  be  assumed  to  be  mathematically  zero  since  this 
violates  the  criterion  of  local  thermodynamic  equilibrium.  The  key 
conclusion  of  the  manuscript  is  that  the  chemical  potential,  /xs,  of 
electrically  neutral  species,  s,  in  a  membrane  need  not  be  mathe¬ 
matically  bounded  by  the  values  at  the  electrodes.  Specifically,  it 
was  shown  that  if  the  ionic  and  the  electronic  currents  through 
the  membrane  are  in  the  same  direction,  the  fis  in  the  membrane 
may  lie  outside  the  range  covered  by  the  values  at  the  electrodes. 
Under  some  conditions,  the  thermodynamic  stability  range  of  the 
membrane  may  be  exceeded  thus  causing  degradation.  It  was  also 
shown  that  during  charging  of  an  electrochemical  device,  ionic  and 
electronic  currents  through  the  electrolyte  are  in  the  same  direc¬ 
tion.  Under  some  conditions  of  charging,  degradation  of  the  device 
may  occur.  In  lithium-ion  cells,  for  example  during  charging,  the 
chemical  potential  of  lithium,  /xLj(x),  in  the  electrolyte  may  exceed 
that  corresponding  to  pure  lithium,  /x°.  This  will  lead  to  the  precip- 
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itation  of  lithium  and/or  reaction  with  the  electrolyte,  depending 
upon  the  operating  conditions  and  the  magnitudes  of  the  various 
transport  resistances.  Capacity  fading  according  to  this  model  thus 
will  primarily  occur  during  charging.  In  battery  packs  containing 
several  cells  connected  in  series,  a  similar  phenomenon  can  occur 
in  an  unbalanced  cell  even  during  discharge.  In  unbalanced  cells  in 
battery  packs,  SEI  layers  may  form  at  both  anode/electrolyte  and 
cathode/electrolyte  interfaces.  If  the  cathode  is  intrinsically  unsta¬ 
ble  in  contact  with  electrolyte,  then  it  will  mean  that  the  SEI  layer  at 
the  cathode  will  be  thicker  in  unbalanced  cells  compared  to  those 
operating  normally.  Thus,  in  unbalanced  cells  in  a  battery,  capacity 
fade  may  occur  during  both  charge  and  discharge.  The  present  work 
also  suggests  that  a  bi-layer  separator  with  one  layer  exhibiting 
high  electronic  conductivity  (placed  in  contact  with  the  cathode) 
may  help  lower  the  /xLi(x)  in  the  electrolyte,  and  thereby  suppress 
degradation.  This  is  a  radically  different  approach  from  the  current 
practice.  It  is  suggested  that  by  suitably  tailoring  transport  prop¬ 
erties  of  the  electrolyte/separator,  the  SEI  layer  and  the  interfaces, 
incidence  of  degradation  may  be  avoided  or  at  least  suppressed. 
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Appendix  A.  Coupled  motion  of  spheres  in  a  fluid 


Fig.  2(a)  shows  two  spheres  of  radii  r\  and  r2  immersed  in  a 
fluid  of  viscosity  ??,  joined  by  a  mass-less,  infinitely  thin  string  going 
over  a  frictionless  pulley  (pulley  outside  the  fluid).  Fig.  2(b)  shows 
two  spheres  immersed  in  a  viscous  fluid  connected  by  a  string,  but 
without  a  pulley.  Apparent  masses  of  the  spheres,  corrected  for 
buoyancy  are  mi  =  {47ir\/3){pM  -  Pfi)g  and  m2  =  (47rr|/3)(pM  - 
Pfi)g  where  pM  and  pp  are  specific  gravities  of  the  spheres  and 
the  fluid,  respectively.  Here  we  have  chosen  mi  >m2.  Fig.  2  shows 
ri  >r2.  If  the  specific  gravities  of  the  two  spheres  are  different,  for 
example  if  they  are  made  of  two  different  materials,  then  it  is  pos¬ 
sible  to  have  mi  >  m2  but  ^  may  be  larger  or  smaller  than  r2.  Note, 
ri  and  r2  are  contained  in  the  transport  (kinetic)  coefficients;  while 
mi  and  m2  are  contained  in  the  thermodynamic  driving  forces. 
In  Fig.  2(a)  the  heavier  sphere  moves  down  (lowers  its  potential 
energy)  while  the  lighter  one  rises  (increases  its  potential  energy). 
There  is  a  net  decrease  in  potential  energy,  which  is  dissipated  as 
heat  lost  to  the  reservoir  (and  eventually  to  the  surroundings).  The 
velocities  of  the  two  spheres  in  steady  state  (terminal  velocities 
achieved)  are  given  by 


V\  =  -v  = 


1 


67H7(ri  +r2) 


(-mig)- 


1 


67H7(ri  +r2) 


{~m2  g) 


=  -TnV/xi  -Li2V/x2  =LnXi  +Li2X2 


(Al) 


The  coupling  tension  in  the  string  is 

(m) 

We  can  also  write  for  the  forces 

Xi  = -Vam  = -Vjt  2-Vl=X2-Vl=X2-Rv  (A5) 

Note  that  LnL22  -l\2  =  0  and  notLnL22  - l\2  >  0,  and  represents 
fully  coupled  processes  (naturally,  as  the  spheres  are  connected  by 
a  string).  In  this  case,  the  forces  are  independent  but  the  fluxes  are 
dependent.  This  means  we  may  choose  Xi  and  X2  independently 
(by  independently  choosing  the  masses  of  the  two  spheres),  but 
the  fluxes  (dependent)  are  given  by  v\  =  -v  =  -v2. 

The  dissipation  function  is  [10] 

-  m2f  >  0 

(A6) 


where  d/S/dt  is  the  entropy  generation  rate. 

Now  consider  a  case  in  which  there  is  no  pulley  and  both 
spheres  (still  joined  by  a  string)  are  falling  under  gravity  in  the 
fluid  (Fig.  2(b)).  Now  the  velocities  are 

v\  —  -v  =  - — — * - -(-mig)  +  - — — * - -(-m2g) 

67H7(ri  +r2)  6jtr]{r i  +  r2y 

=  -f-nV/Xi -T12V/x2  =L11X1  +T12X2  (A7) 

and 

=  -^21  V/Xi  -  I22V/x2  =  L2iX!  +L22X2  (A8) 


Again  the  forces  are  independent  but  the  velocities  (fluxes)  are 
dependent.  Now 

<*>) 

Note  the  coupling  coefficients  are  positive,  that  is 


^12  =  ^21 


1 

67H7(ri  +r2) 


>  0 


(A10) 


Again  LUL22  - L\2  =  0,  consistent  with  fully  coupled  processes. 
Coupling  tension  in  the  string  is 


(A,i> 

We  can  also  write  for  the  forces 

X,  =  -V/i!  =  WH2-Vl=-X2-Vl=  — X2  -  Rv  (A12) 

The  dissipation  function  is 


and 

"2  -  "  —  6»,,(rl+r2)<— "'S’  +  55S7W-'"!S) 
=  -^2iV/Xi  -I22V/x2  =L2iX!  +L22X2 


^  =  T^-  =  YjiX,  =  ~v\  VAM  -  ViVn2 


(A2) 


gz  9 

6  (mi  +  m2y  >  0 


67r(ri  +r2) 


(A13) 


The  above  are  velocities,  so  the  units  are  ms-1.  The  coupling 
coefficients  are  negative,  that  is 


^12  =  f-21 


1 

67H7(ri  +r2) 


<  0 


(A3) 


This  case  corresponds  to  Ln  =L22  =  -I12  =  -L21  =1  =  1/ 
(67H7(ri  +r2))>0. 


In  this  case,  both  spheres  move  down  (both  reduce 
their  potential  energies).  The  coupling  tension  is  zero  if 
(r2/(r  i  +r2))mi  =(n/(ri  +r2))m2. 

The  signs  of  the  forces  are  the  same  in  both  cases,  but  the  signs 
of  L12  =  L2\  are  different.  Thus  in  the  two  cases  the  signs  of  the  terms 
Li2X2  and  L2iXi  are  different.  In  one  case,  in  which  Ii2X2  >  0  and 
L2iXi  >  0,  an  increase  in  energy  of  one  sphere  occurs  as  a  result  of 
coupling. 
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Appendix  B.  Transport  through  predominantly  ionic 
conductors 

In  liquid  electrolytes  containing  ionic  salts  dissolved  in  polar 
liquids,  there  is  at  least  one  mobile  cation,  and  at  least  one  mobile 
anion.  For  an  electrolyte  MmZn,  a  salt  of  metal  M  and  a  nonmetal  Z, 
local  equilibrium  gives 

MmZn(r)  ^  mMn+(r)  +  nZm-(r)  (B1 ) 


Thus,  if  electrons  and  Zr4+  (or  Zr)  exhibit  nonzero  conductiv¬ 
ities  (diffusivities),  local  equilibrium  criterion  can  be  satisfied.  In 
general,  we  expect  sufficient  (nonzero)  mobilities  of  the  various 
species  to  satisfy  local  equilibrium. 

We  now  consider  local  equilibrium  at  r  between  neutral  species 
s,  the  corresponding  constituent  mobile  ions  i  of  valence  z2,  and 
electrons,  e,  namely 

s(r)  i(r)  +  Zj-e(r)  (B12) 


and 

=  m/xMn+(f)  +  n^Zm-(r)  =  mjlMn+{f)  +  nfiZm-{r)  (B2) 

where  A;(? )  =  /x2(r)  +  z2e0(r)  is  the  local  electrochemical  potential 
of  i,  Zi  is  the  valence  of  i,  e  is  the  electronic  charge,  and  @(r )  is  the 
local  electrostatic  potential.  Local  equilibrium  also  implies 

M(r)  <£>  Mn+(r)  +  ne(r)  (B3) 

which  gives 

mtf)  =  A  M"+  (?)  +  nAe(?)  (B4) 

and 

Z(r)  +  me(r)^Zm-(r)  (B5) 

which  gives 

Mz(?)  =  A zm-  (?)  -  mAe(?)  (B6) 


where  Ae(?)  =  /xe(r)  -  e<2>(?)  is  the  electrochemical  potential  of 
electrons.  In  salt  solutions  in  polar  liquids,  free  electron  concentra¬ 
tion  is  very  small  (which  could  be  present  as  solvated  electrons). 
Thus,  negligible  electron  current  transports  through  the  solution. 
An  important  point  is  that  chemical  and  electrochemical  potentials 
of  electrons  are  still  well-defined  quantities  in  salt  solutions,  some 
concentration  of  free  or  solvated  electrons  always  exists,  and  some 
electron  transport  must  also  occur  when  nonzero  electrochemical 
potential  gradient  of  electrons  exists.  Finally,  for  local  equilibrium 
to  exist,  differential  forms  of  the  above  equations  must  also  be  valid, 
that  is, 

<$MMmz„(?)  =  m8jLMn+{r)  +  n8/nZm-(r)  =  m8jlMn+(r)  +  mSAz™-(?) 

=  +  nSfizir)  (B7) 

and 

<$Mm(?)  =  <W(?)  +  n<$Ae(?)  (B8) 

where  8 ft  or  8jl  denote  differential  variations  in  chemical  potential 
fi  or  electrochemical  potential  A-  For  a  line  compound  (a  very  nar¬ 
row  range  of  stoichiometry),  8fiMmzn(j)  =  0-  Even  then,  individual 
chemical  potentials  of  neutral  metal  and  neutral  nonmetal  species 
may  exhibit  large  variations. 

For  differentials  to  satisfy  local  equilibrium,  relevant  species 
must  exhibit  some  transport;  that  is  a  sufficient  number  of  species 
comprising  the  system  must  exhibit  nonzero  conductivities  or  dif¬ 
fusivities.  In  typical  solid  electrolytes,  there  often  exists  only  one 
substantially  mobile  ion.  Even  so  the  assumption  of  local  equilib¬ 
rium  implies  some  mobility  for  a  sufficient  number  of  species  to 
satisfy  local  equilibrium.  For  example,  consider  yttria-stabilized 
zirconia  (YSZ)  which  is  a  predominantly  oxygen  ion  conductor. 
Local  equilibrium  for  neutral  Zr  is  given  by 

Zr(r)  Zr4+(r)  +  4e(r)  (B9) 

For  local  equilibrium  to  be  satisfied,  we  must  have 
Mzr(r)  =  /U,Zr 4+(r)  +  4  fie{r)  =  AZr4+(t)  +  4  Ae(?)  (BIO) 

and 

8n  zr(?)  =  <$MZr4+(?)  +  4  8  MO  =  SjlZr4+(r)  +  4Sjle(f)  (B1 1 ) 


with 

r  ai(f)(ae(r))z‘ 

K=  as(f)  (B13) 

as  the  equilibrium  constant  and  a's  are  local  thermodynamic  activ¬ 
ities  of  species  in  Eq.  (B12). 

The  valence  z2-  can  be  positive  (cation)  or  negative  (anion).  Local 
equilibrium  at  r  gives 

Ms(r)  =  Hi(j)  +  Ziflelr)  =  Ai(?)  +  z,Ae(?)  (B14) 


where  /xs(r),  /xe(?)  and  /x2(r)  are  respectively  the  chemical  poten¬ 
tials  of  neutral  species  (s),  electrons  (e),  and  ions  (i)  at  r;  and  Ai(?) 
and  Ae(?)  are  respectively  the  electrochemical  potentials  of  ions 
and  electrons. 

Experimentally  measurable  electric  potential  is  given  by  [33,34] 


<p(r)  = 


Ae(r)  =  _Mr)  + 

e  e 


(B15) 


The  chemical  potential  of  neutral  species  s  is  thus  also  given  by 


/ts(r)  =  A i(?)  +  z,Ae(?)  =  Mr)  ~  Zie<p{r) 


(B16) 


The  above  refers  to  local  equilibrium,  which  means 

A Ais(r,  r0)=/u,s(?)  -  fis(r0 )=M?)  -  Mfo)  +  z,{Ae(?)  -  Ae(?o)}  (B17) 


or 


AMs(r,  r0)  =  Ms(?)  -  Mh)  =  Mr)  -  Ai(?0)  —  z,e{^(r)  -  y>(f0))(B18) 

as  the  difference  in  chemical  potentials  of  s  between  r  and  f0. 

Eqs.  (B17)  and  (B18)  are  the  result  of  local  equilibrium.  But  we 
also  have  from  the  Onsager  equations 

A  Mr,  f0)  =  M?)  -  Mfo) 


RskJk  ■  dr 


where  k  =  s,  i,  e 


(B19) 


which  means 


A  Mr, f0)  =  MO -Mfo) 


=  -  /  RssJs  df-f  RsJi  dr  -  f  RseJe  ■  dr 

«7  To  J  to  J  To 


(B20) 


Strictly,  Eqs.  (B19)  and  (B20)  should  include  all  other  species  in 
the  system  also.  The  neglect  of  the  other  species  is  an  approxima¬ 
tion.  The  ionic  and  electronic  fluxes  are  respectively  j2(r)  and  Je(r) 
in  #m-2  s-1  and  can  be  written  in  terms  of  the  respective  current 
densities  as 

;,(?)  =  (B2i ) 

z2e 

and 

]e(f)= (B22) 


Substituting  for  fluxes  in  terms  of  current  densities  in  Eq.  (B20) 
gives 


A/xs(r,  r0)  =  /xs(r)-/xs(r0) 
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—  —  f  RssJs  •  dr  —  —  f  Rsil  •  dr2  +  —  f  Rsele  •  dr 
Jfo  Zl‘e  Jr0  e  Jr0 


(B23) 


We  will  now  assume  that  the  concentration  of  the  neutral 
species,  s,  within  the  material  is  negligible.  This  means  a  very  high 
equilibrium  constant  for  the  reaction  soi  +  z2e.  Thus,  there  is  hardly 
any  concentration  of  neutral  s;  that  is  Cs(r)  ~  0  and  thus  hardly  any 
transport  of  neutral  s  (as  neutral  s)  occurs  through  the  material; 
that  is  Js(j)  =  0.  If  the  material  is  a  predominantly  ionic  conductor, 
it  also  means  free  electron  concentration  is  very  small.  However,  we 
still  cannot  set  Ie  identically  zero  to  ensure  that  local  equilibrium  is 
satisfied.  Thus,  Eq.  (B23)  reduces  to 


A/xs(r,  r0)  =  Ms(r)-/xs(r0) 


1 

z^e 


RSih '  dr  + 


;  /*• 

Jr0 


le  •  dr 


(B24) 


is  the  ionic  area  specific  resistance  between  x'0  and  x,  and 

r'(x)  =  (x~x«)  (B31 ) 

&e 

is  the  electronic  area  specific  resistance  between  x'0  and  x. 

Thus, 


Ms(x)  =  Hs{x'o)-zier'i{x)Ii  +  zier'e(x)Ie 


=  Ms(Xo)  -  zfe  [r-(x)/,-  -  r'(x)/eJ  (B32) 

which  is  in  a  form  similar  to  the  author’s  previous  derivation  [13]. 

Comparison  of  Eqs.  (B26)  and  (B29)  gives  relations  between  the 
Onsager  coefficients  and  conductivities  as 


Rsi  — 
and 
Rse  — 


(B33) 


(B34) 


The  Rss  is  finite,  nonzero,  positive.  However,  by  making  the  usual 
measurements  such  as  currents  and  voltages,  we  cannot  determine 
Rss,  the  parameter  related  to  the  transport  of  neutral  s.  It  can  in  prin¬ 
ciple  be  determined  if  in  addition  to  the  measurements  of  ionic  and 
electronic  currents,  one  also  accurately  measures  loss/gain  of  s  at 
the  two  electrodes  (reservoirs)  which  would  require  the  experi¬ 
ments  to  be  conducted  for  a  sufficiently  long  time.  Tubandt  used 
this  approach  some  90  years  ago  to  study  transport  through  silver 
halides  [35]. 

We  now  consider  one  dimensional  transport  and  assume  that 
the  transport  coefficients  are  independent  of  position.  Referring  to 
Fig.  1,  we  will  initially  choose  reference  point  x'0  within  the  mem¬ 
brane.  Then,  Eq.  (B24)  becomes 


/x  ,  ,  x  *  ,  ..  RrilAx-X')  Rsele(x-x'n) 

fls(x)  -  fls(x'o)  =  A fls(X,  X' )  =  -  51  0>  +  - o_i  (B25) 


z,e 


or 


/x  ,  ,  N  RsMx-X')  Rsele(x-x') 
Ms(x)  =  fisix'o)  ~  51  _  _  °J  +  -eeK  oJ 


z,e 


(B26) 


These  are  the  Onsager  coupling  coefficients  which  relate  forces 
(treated  as  dependent  variables)  to  fluxes  (treated  as  independent 
variables).  Note  also  RSI  /  ooandRse  +  oo  (finite),  indicating  a  cou¬ 
pling  of  fluxes. 

Appendix  C.  Abrupt  changes  across  interfaces 

Let  us  consider  points  just  inside  the  membrane  adjacent  to  the 
electrodes  (reservoirs)  given  by  x  =  x0  +  8Xo  and  x  =  Xf  -  8Xf,  as  in 
Fig.  1.  Note  that  8Xo ,  8Xf  «  jx/  -  x0 1 .  The  chemical  potentials  of  s  in 
the  membrane  close  to  the  interfaces  are  given  by  [using  Eq.  (B32)] 

Ms(x0  +  SXo )  =  Ms(x0)  -  z,e  [rf°Ii  -  r^°/eJ  =  Ms(x+ )  (Cl ) 

and 

fis{Xf  -8Xf)  =  fjLs(xf)  +  zte  [rffli  -  r*fle\  =  fis{x-f)  (C2) 

where  xj  =  x0  +  is  slightly  greater  than  x0  and  Xj  =  Xf  -  8Xf  is 

slightly  less  than  Xf.  The  corresponding  (abrupt)  changes  in  fis 
across  the  interfaces  are 


Since  we  are  considering  a  one  dimensional  steady  state  prob¬ 
lem,  If  and  Ie  are  position  independent.  Units  of  RS1  and  Rse  are 
erg  cm  s  or  Jms  (with  chemical  potential  given  per  species— per 
ion,  per  electron,  per  atom,  per  molecule). 

No  assumptions  concerning  the  relationships  between  current 
densities  and  thermodynamic  potential  gradients  (such  as  electro¬ 
chemical  potential  gradients)  have  been  made  so  far.  Let  us  now 
assume  that  the  flux  of  a  charged  species  occurs  down  its  electro¬ 
chemical  potential  gradient.  That  is, 

/,(?)  = and  ?e(?)-£v£e(?)  =  -<reVtfr)  (B27) 

where  a ,•  and  cre  are  respectively  the  ionic  and  the  electronic  con¬ 
ductivities.  Eq.  (B27)  represents  an  important  assumption  made  in 
virtually  all  transport  studies  of  electrically  charged  species  [36]. 
Substitution  in  Eq.  (B18)  for  a  one  dimensional  problem  gives 

/*s(x)-Ms(x'0)=_zg  zje 

v  _  v'  nr.  1  v  ' 


or 


£ts(x)  -  Ats(x'o)  = 


Zje(x-x'0)r  t 
Gi 


z,e(x  -  Xp) 

Oe 


Ie 


(B29) 


Now 

r-(x)  =  (B30) 

1  O’,- 


A nxs°  =  Zj-e  |  rf°li  -  rxe°Ie\ 

(C3) 

and 

A  Hxsf  =  Zfi  (rff  Ij  -  if  Ie] 

(C4) 

Note  in  general,  =  z{e  ^rx°/2-  -  rx°/eJ  j=  0  and  A/xx/  = 

z2-e  [rffli  -  r%f  I e]  ^0.  With  fis{Xf)  j=  fis{x0)  (experimentally 
imposed  conditions),  the  only  situation  for  which  A/xXo  =  0, 
corresponds  to  both  rXo  =  0  and  rx°  =0;  and  the  only  situation  for 
which  =  0,  corresponds  to  both  rff  =  0  and  r*f  =  0.  These 
are  very  unlikely  situations  since  they  imply  perfect  reversibility 
for  the  transfer  of  both  ions  and  electrons  across  interfaces.  Thus, 
A/xXo  and  A/xx/,  which  represent  sharp  changes  in  /xs  across  elec¬ 
trode  ( reservoir )/electrolyte  interfaces  are  nonzero.  Let  us  assume 
that  the  electrodes  are  perfectly  (electrochemically)  reversible,  a 
terminology  commonly  used  in  electrochemical  literature.  This 
means  r?°  =  0  and  rff  =  0.  Even  in  such  a  case  note  that 

\An*>\  =  \zfiIel?\±0  (C5) 

and 

\Afixsf\  =  |z2e/erx/|  /0  (C6) 

The  preceding  shows  that  a  perfectly  reversible  electrode  in 
an  electrochemical  sense  does  not  imply  the  equilibration  of 
fis,  the  chemical  potential  of  neutral  species,  across  an  elec¬ 
trode/electrolyte  interface  (when  /xs(Xf)  /  /xs(x0)),  and  that  there 
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will  in  general  be  abrupt  changes  in  /xs  across  interfaces.  Per¬ 
fectly  reversible  electrode  in  an  electrochemical  sense  only  means 
the  electrochemical  potential,  jli}  of  the  mobile  ion  equilibrates 
across  the  interface  even  at  finite,  nonzero  current  densities.  There 
is  considerable  literature  in  electrochemistry  which  incorrectly 
assumes  the  equilibration  of  the  chemical  potential  of  neutral 
species,  /zs,  across  electrode/electrolyte  interfaces  when  electrodes 
are  reversible  in  an  electrochemical  sense. 
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